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Abstract 

One interesting method of constraining the dense matter Equations of State is to 
measure the advancement of the periastron of the orbit of a binary radio pulsar (when it 
belongs to a double neutron star system). There is a great deal of interest on applicability 
of this procedure to the double pulsar system PSR J0737-3039 (A/B). Although the above 
method can be applied to PSR A in future within some limitations, for PSR B this method 
can not be applied. On the other hand, the study of genesis of PSR B might be useful in 
this connection and its low mass might be an indication that it could be a strange star. 
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1 Introduction 



It was argued by Witten [1] that strange quark matter, composed of u, d and s quarks, is more 
bound than nuclei at zero pressure and is therefore the true ground state of matter. Ever since 
the advent of this hypothesis the existence and the observational characteristics of the strange 
quark stars [2, 3] are being debated. 

Unfortunately, strange stars and neutron stars have a number of almost similar physical 
properties like the range of masses, stable rotation periods, cooling history and so on, but 
strange stars are more compact than neutron stars. Because of this, quite a few compact 
objects were speculated to be strange stars but there has been no conclusive evidence for the 
case or against it. For example, RX J1856. 5-3754 has been the subject of such a recent debate 
because of certain peculiarities of its thermal X-ray properties [4, 5]. 

Understandably the most convincing proof for a strange star would come from a direct 
indication of the equation of state and we believe that the double neutron systems (DNS) 
might provide us with that answer. Precise measurements of a number of parameters of DNS 
systems allow for the study of a variety of phenomena like the effects of general relativity, the 
emission of gravitational waves or the intricacies of the binary evolution. In particular, in the 
case of PSR J0737— 3039, the only double pulsar system [6, 7] detected so far, it is possible to 
measure many more properties like the mass ratio of the two neutron stars directly obtainable 
from pulsar timing. The two pulsars are named as J0737— 3039A and J0737— 3039B. The values 
of spin period and magnetic field of these two are P = 22.7 ms, B = 6.4 x 10 9 G for pulsar A 
and P = 2.77 s and B = 1.6 x 10 12 G for pulsar B. Apparently A is a typical recycled pulsar 
and B is the one born second. The periods Pa,b and their time derivatives provide upper limits 
for the lifetimes tA ~ 210 Myr and ts ~ 50 Myr. The estimated separation (sum of the semi 
major axes) was R ~ 10 12 cm and the eccentricity was e ~ 0.11 when pulsar B was born. 
The present values of separation and eccentricity are 8.8 x 10 11 cm and 0.0877775(9). As the 
values of R and e does not change considerably with time, the analysis of the origin of PSR 
J0737— 3039B is insensitive to the exact age of the pulsar. 

Formation models of PSR J0737— 3039 were considered by various authors [8, 9, 10]. Dewi 
& van den Heuvel [8] proposed that J0737— 3039 has originated from a close helium star plus 
neutron star (HeSNS) binary in which at the onset of the evolution the helium star had a mass 
in the range 4.0 — 6.5 M and an orbital period in the range 0.1 — 0.2 day. In their model, a 
kick velocity in the range 70 — 230 km /sec must have been imparted to the second neutron 
star at its birth. Podsiadlowski et al. [9] proposed an electron capture supernova scenario. 
Using the observational indications that the binary system will end up within 50 pc from the 
galactic plane with 0.088 < e > 0.14 and with a transverse velocity less than 30 km/s, provided 
50 Myrs in the past the progenitor system had a circular orbit, Piran & Shaviv [10] estimated 
that it is most likely that the progenitor of PSR J0737-3039B had a mass of 1.55 ± 0.2 M 
with a low kick velocity 30 km/s. The progenitor may well have been a NeOMg star. Willems 
et al. [11] used population synthesis method to estimate the most likely radial velocity of the 
binary. They favored the standard formation scenario with reasonably high kick velocities (70- 



2 



180 km/s) and pre-SN mass m^i of at least 2 M which fits with a Helium star progenitor, but 
acknowledged that for very small transverse velocities (~ 10 km/s) lower mass progenitors are 
allowed or even favored. 



2 The nature of the members of PSR J0737-3039 

The system is rich in observational phenomena, including a short radio eclipse of A by B and 
orbital modulation of the radio flux of B affected by A. Much radio pulse data from both 
pulsars has already been accumulated. High energy observations could also be important to 
these studies since both stars also emit X-rays and they are now being detected recently. 



2.1 radio observations 

The best answer to our question regarding the equation of state would come from the measure- 
ment of the moment of inertia of individual stars in this system. The moment of inertia can 
be determined by measuring the extra advancement of the periastron of the orbit above the 
standard post-Newtonian one due to the spin-orbit coupling [12] using timing analysis of radio 
data. Periastron advance (which is the observable quantity) for a DNS system is given by : 
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where the first term inside the square bracket represents 1PN term, the second term term inside 
the third bracket represents 2PN term, and the third term inside the third bracket represents 
spin-orbit coupling (SO). The parameters used in the above expression are defined as : 
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where G is the gravitational constant, I a is the moment of inertia of a th body (a = 1, 2), 
uj a = 2n /P spin is its angular velocity of rotation (here a + 1 means modulo 2, i.e. 2+1=1) 



3 



and n = 2tt / P or b- %i = mi/M, x 2 = m^jM and M = mi + 777,2, k is the orbital angular 
momentum vector, s a is the spin vector and Ko is the line of sight vector. Using eqn 1, the 
moment of inertia can be estimated from the observed value of Co provided we know P sp i n , e, 
Porb, m i and m 2 (which can be estimated from timing analysis of radio data) of the binary 
components which is the case for J0737-3039 A and B. 

On the other hand, the moment of inertia of a star of known mass can be theoretically 
calculated using a particular EoS. Bejger et al. [13] calculated the moment of inertia for both 
PSR J0737— 3039 A and B using 25 EoS of dense matter including an EoS for strange quark 
matter. Bejger et al. [13] found I a to lie between 1.894 x 10 45 g cm 2 and 0.726 x 10 45 g cm 2 . They 
also checked that the moment of inertia estimated using Hartle's approximations (as described 
in Kalogera & Psaltis [14]) differs from the exact calculation with rotating axisymmetric stellar 
model by not more than 1%. We re-performed the analysis and found I a to lie between 
0.67 - 0.74 x 10 45 g cm 2 . It is worthwhile to mention that the stiffest EoS (MFT17) chosen 
by Bejger et al. strongly overestimate and the softest EoS (BPAL12) strongly underestimate 
nuclear matter incompressibility. So from the view point of a theoretician, the value of I a is 
likely to be somewhere in the middle region of figure 1 of Bejger et al. [13]. 

But the point to note here is that, in the expression of Co, the moments of inertia of both 
members of a DNS system are present (eqn 1) which means that we have one observable and two 
unknowns which is impossible to solve. But we have checked that the f3 sa term is significant only 
for P sp in < 100 ms, so the term for J0737— 3039B can be neglected and only the contribution 
from J0737— 3039A remains. So determining the moment of inertia of PSR J0737— 39B is 
impossible where measuring the same for PSR J0737— 39A is possible by neglecting the spin- 
orbit coupling term due to PSR J0737— 39B. This implies that this method of understanding 
the form of the matter building the star from pulsar timing analysis can be used for PSR 
J0737— 39A, but not for PSR J0737— 39B. So we can not test the prediction of I B as claimed 
by Bejger et al. [13]. 

Even for PSR A, to measure the moment of inertia within 10% of accuracy, we need to 
improve the accuracy in the knowledge of binary parameters to a significant extent. This 
possibility was studied recently by Iorio [15] who concludes that "the possibility of reaching in 
a near future the required accuracy to effectively constrain I a to 10% level should be, perhaps, 
considered with a certain skepticism" . 

2.2 genesis of PSR J0737-39B 

One can consider an estimate of minimum progenitor mass for PSR J0737— 39B, keeping in 
mind the recent comment by Stairs et al. [16], that an important point usually neglected is the 
contribution from gravitational binding energy which is defined as 

E grav = M P -M G (6) 
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where Mq is the gravitational mass and Mp is the proper mass. The increase in baryonic 
energy is given by 

r* AE(r)n(r)dr 

Jo N _ 2GM(r) l V2 



bar 



where AE is the increase in energy per baryon due to the conversion from ONeMg matter to 
strange quark matter, n is the number density, G is the gravitational constant, M and R are 
the mass and radius of the star. So the mass of the progenitor will be 

Mp r0 g en itor = M pu [ sar + E tot + Am (8) 

where 

Etot = E grav + Eb ar (9) 

and Am is the mass loss during the conversion process. PSR J0737— 3039B, which has M G = 1.249M , 
corresponds to Mp = 1.5187 in our SS model [17] and we found E grav = 0.2697M Q and 
E bar = 0.0708M© giving E tot = 0.3405M Q . So M progenitor = 1.249 + 0.3405 + Am = 
1.5895 + Am. This value is close to that estimated by Piran and Shaviv [10] i.e. 1.55 ± 0.2 M Q 
with a low kick velocity 30 km/s. On the other hand, if J0737— 3039 is a neutron star, 
we get using APR EoS [18], M P = 1.4357 giving E grav = 0.1867 and E bar = giving 
M progen itor = 1-249 + 0.1867 + Am = 1.4357 + Am which is also not far from Piran and 
Shaviv's. Remember, detailed simulation is needed to know the exact value of Am which we 
hope will be done in the future. So it is very difficult to identify whether J0737 — 3039 B is a 
strange star or a neutron star from these arguments. 

Another interesting outcome is that, in our model [17], a strange star with Mq = 1.249 M & 
has 1.9058 x 10 57 baryons while a neutron star with M G = 1.249 M has 1.6390 x 10 57 while a 
normal NeOMg star of mass 1.59 M & has 1.90 60 x 10 57 number of baryons. We have assumed, 
of course, that the gravitational mass is the same as its baryonic mass for a normal star - as 
the general relativistic effects are very small. This estimate shows that if J0737 — 3039 B is 
a strange star then the baryon loss is negligible (0.01%) implying a small Am. On the other 
hand, if J0737 — 3039 B is a neutron star then the baryon loss is larger (14%) implying larger 
values for Am. Small baryon loss (as obtained in the first case) supports Piran and Shaviv's 
model with small kick. So this study would suggest PSR B to be a strange star. 



3 Conclusions and summary 

We have discussed how the double pulsar system can help us to constrain the dense matter 
EsoS. Radio observations (timing analysis) may help us to constrain the EoS of PSR A in future 
and may allow us to deduce whether it is a strange star or a neutron star. 

For PSR B, as this is not possible in the present situation, we can use indirect methods, e.g. 
modeling its formation scenario from the present day observed parameters. The small mass of 
the star might be an indication that it is a strange star, since in the formation of such a star, 
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in addition to the gravitational binding energy loss, mass is also lost in the form of baryonic 
energy. Although this is suggestive, more modeling is required to support such a conclusion. 
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